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^^ Abstract. Even at weak to intermediate coupling, the Hubbard model poses a 

T-H formidable challenge. In two dimensions in particular, standard methods such as 

the Random Phase Approximation are no longer valid since they predict a finite 

' I ' temperature antiferromagnetic phase transition prohibited by the Mermin- Wagner 

D theorem. The Two-Particle-Self-Consistent (TPSC) approach satisfies that theorem 

as well as particle conservation, the Pauli principle, the local moment and local 

charge sum rules. The self-energy formula does not assume a Migdal theorem. There 

is consistency between one- and two-particle quantities. Internal accuracy checks 

5^ allow one to test the limits of validity of TPSC. Here I present a pedagogical review 

ri of TPSC along with a short summary of existing results and two case studies: a) 

I the opening of a pseudogap in two dimensions when the correlation length is larger 

'■^ than the thermal de Broglie wavelength, and b) the conditions for the appearance 

^ of d-wave superconductivity in the two-dimensional Hubbard model. 
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1.1 Introduction 

Very few models can describe complex behavior observed in nature with an 
economy of parameters. The Hubbard model is in this category. It has be- 
come the cornerstone of correlated electron physics. On the down side, it is 
extremely difficult to solve. While it was proposed in 1963 [1-3], the only exact 
results that we know are in one dimension [4] and in infinite dimension [5]. A 
variety of approximate approaches to solve this model exist, as can be checked 
from the table of contents of this volume. The Two-Particle-Self-Consistent 
(TPSC) [6-8] approach that is described in the present Chapter is in the 
category of non-perturbative semi-analytical approaches. By semi-analytical, 
I mean that while it is possible to find many analytical results, numerical 
integrations are necessary in the end to obtain quantitative results. 

Why should you bother to learn yet another approach? Because in its 
known regime of applicability it is extremely reliable, as can be judged by 
benchmark Quantum Monte Carlo (QMC) calculations. Because it satisfies 
a number of exact results that control the quality of the approximation and 
make it physically appealing. And because it gives physical insight into many 
questions related to the two-dimensional Hubbard model relevant for the high- 
temperature superconductors and many other materials. As a case study, I 
discuss in this Chapter the physics of pseudogap induced by precursors to 
long-range order. We will see that this describes the physics of the pseudo- 
gap in electron-doped high-temperature superconductors where predictions of 
TPSC have been verified experimentally. Pseudogap phenomena include the 
appearance of a minimum in the single particle spectral weight and density 
of state at the Fermi level. I will be more precise later. 

This Chapter offers to the reader a simple pedagogical introduction to this 
approach along with the case studies mentioned above and a guide to various 
other problems that have been, or have not yet been, solved with TPSC. 

I assume familiarity with the basics of many-body theory, i.e. the canonical 
formalism, second quantization, many-body Green's functions, response func- 
tions and with the Matsubara formalism for finite temperature calculations. 
Knowledge of functional derivative approaches would be useful for some of the 
more advanced topics, but it is not essential to learn the important results. 

Before you read on, you might be interested to know a little more about 
the method to decide whether it is worth the effort. TPSC is designed to study 
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the one-band Hubbard model 

i? = - ^ U^jcla^ja + UY^ n,^n,^ (1.1) 

ij(T i 

where the operator Ci<j destroys an electron of spin a at site i. Its adjoint c\^ 
creates an electron and the number operator is defined by ni„ — cl^Ci^. The 
symmetric hopping matrix ti_j determines the band structure, which here can 
be arbitrary. The screened Coulomb interaction is represented by the energy 
cost U of double occupation. It is also possible to generalize to cases where 
near-neighbor interactions are included. We work in units where fc^ = 1, ft = 1 
and the lattice spacing is also unity, a = 1. In all numerical calculations, we 
take as unit of energy the nearest-neighbor hopping t — I. 

One of the first concepts that is discussed with the Hubbard model is that 
of the Mott transition [9]. When dimension is larger than unity, at "strong 
coupling", large U/t, the states are localized, but at "weak coupling", small 
U/t, the states are delocalized. The Mott transition is quite subtle and has 
been the subject of many papers. It is discussed in the Chapters of M. Jarrell, 
M. Potthoff, D. Senechal and others in this volume. 

TPSC is valid from weak to intermediate coupling. Hence, on the negative 
side, it does not describe the Mott transition. Nevertheless, there is a large 
number of physical phenomena that it allows to study. An important one is an- 
tiferromagnctic fluctuations in two- or higher-dimensional lattices. A standard 
Random Phase Approximation (RPA) calculation of the spin susceptibility 
signals a finite temperature phase transition to antifcrromagnetic long-range 
order. This is prohibited by the Mcrmin- Wagner theorem [10,11] that states 
that in two dimensions you cannot break a continuous symmetry at finite tem- 
perature. It is extremely important physically that in two dimensions there is 
a wide range of temperatures where there are huge antiferromagnetic fluctu- 
ations in the paramagnetic state. The standard way to treat fluctuations in 
many-body theory, RPA misses this. As we will see, the RPA also violates the 
Pauli principle in an important way. The composite operator method (COM), 
described in this volume by Avella and Mancini (see Chap. ??), is another 
approach that satisfies the Mermin- Wagner theorem and the Pauli princi- 
ple [12-14]. What other approaches satisfy the Mermin- Wagner theorem at 
weak coupling? The Fluctuation Exchange Approximation (FLEX) [15,16], 
and the self-consistent renormalized theory of Moriya-Lonzarich [17-19]. Each 
has its strengths and weaknesses, as discussed in Refs. [7,20]. Weak coupling 
renormalization group approaches "^ become uncontrolled when the antiferro- 
magnetic fluctuations begin to diverge [21-24]. Other approaches include the 
effective spin-Hamiltonian approach [25]. 

In summary, the advantages and disadvantages of TPSC are as follows. 
Advantages: 

^ See the contribution of C. Honerkamp in this volume. 
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There are no adjustable parameters. 

Several exact results are satisfied: Conservation laws for spin and charge, 

the Mermin- Wagner theorem, the Pauli principle in the form ( n2 \ = (n^) , 

the local moment and local-charge sum rules and the f sum-rule. 
Consistency between one and two-particle properties serves as a guide to 
the domain of validity of the approach. (Double occupancy obtained from 
sum rules on spin and charge equals that obtained from the self-energy 
and the Green function). 

Up to intermediate coupling, TPSC agrees within a few percent with Quan- 
tum Monte Carlo (QMC) calculations. Note that QMC calculations can 
serve as benchmarks since they are exact within statistical accuracy, but 
they are limited in the range of physical parameter accessible because of 
the sign problem. 

We do not need to assume that Migdal's theorem apphcs to be able to 
obtain the self-energy. 

The main successes of TPSC that I will discuss include 

Understanding the physics of the pseudogap induced by precursors of a 
long-range ordered phase in two dimensions. For this understanding, one 
needs a method that satisfies the Mermin- Wagner theorem to create a 
broad temperature range where the antiferromagnetic correlation length 
is larger than the thermal de Broglie wavelength. That method must also 
allow one to compute the self-energy reliably. Only TPSC does both. 
Explaining the pseudogap in electron-doped cupratc superconductors over 
a wide range of dopings. 

Finding estimates of the transition temperature for d-wave superconduc- 
tivity that were found later in agreement with quantum cluster approaches 
such as the Dynamical Cluster Approximation. 

Giving quantitative estimates of the range of temperature where quantum 
critical behavior can affect the physics. 

The drawbacks of this approach, that I explain as we go along, are that 

It works well in two or more dimensions, not in one dimension '* [27]. 

It is not valid at strong coupling, except at very high temperature where 

it recovers the atomic limit [28]. 

It is not valid deep in the renormalized classical regime [6]. 

For models other than the one-band Hubbard model, one usually runs out 

of sum rules and it is in general not possible to find all parameters self- 

consistently. With nearest-neighbor repulsion, it has been possible to find 

a way out as I will discuss below. 
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Modifications have been proposed in zero dimension to use as impurity solver for 
DMFT [26] 
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For detailed comparisons with QMC calculations, discussions of the physics 
and detailed comparisons with other approaches, you can refer to Ref. [7,20]. 
You can read Ref. [29] for a review of the work related to the pseudogap 
and superconductivity up to 2005 including detailed comparisons with Quan- 
tum Cluster approaches in the regime of validity that overlaps with TPSC 
(intermediate coupling) . 

Section 1.2 introduces TPSC in the simplest physically motivated way 
and demonstrates the various results that are exactly satisfied. The following 
section 1.3 presents two case studies: the pseudogap and d-wave superconduc- 
tivity. Many more known results and extensions are summarized in section 1.4. 
The attractive Hubbard model is in the next to last section 1.5. I conclude 
with some open problems in section 1.6. 



1.2 The method 

In the first part of this section, I present TPSC as if we were discussing in 
front of a chalkboard. More formal ways of presenting the results come later. 

1.2.1 Physically motivated approach, spin and charge fluctuations 

As basic physical requirements, we would like our approach to satisfy a) con- 
servation laws, b) the Pauli principle and c) the Mermin-Wagner-Hohenberg- 
Coleman theorem. The standard RPA approach satisfies the first requirement 
but not the other two. Let us see this. With the charge and spin given by 

"i = "i^f + n,4 (1.2) 

S: ^ n,^ (r) - n,^ (r) . (1.3) 

the RPA spin and charge susceptibilities in the one-band Hubbard model are 
given respectively by 

with q a short-hand for both wave vector q and Matsubara frequency and 
where Xo{q) is the Lindhard function that in analytically continued retarded 
form is, for a discrete lattice of TV sites, 

x°^(q,-) = -|E^r^^^^^^^'^- (1-5) 

N ^ LU + 17^ + 61^- Ek+q 

In this expression, assuming periodic boundary conditions. 
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with the sum over j is running over ah neighbors of any of the sites i. The 
chemical potential ^ is chosen so that we have the required density. 

It is known on general grounds [30] that RPA satisfies conservation laws, 
but it is easy to check that for a special case. Since spin and charge are 
conserved, then the equalities X%(.'i= O:'^) — ^^^ X^(fl = Oj"^) ~ for 
w 7^ follow from the corresponding equality for the non-interacting Lindhard 
function x°^(q = 0,a;) = 0. 

To check that RPA violates the Mermin Wagner theorem, it suffices to 
note that if U is larger than Uc = 2/xo'(qmax,w = 0), then the denominator 
1 — ■^Uxoio) of Xspio) can diverge at some wave vector qmax and temperature. 

The violation of the Pauli principle requires a bit more thinking. We derive 
a sum rule that rests on the use of the Pauli principle and check that it is 
violated by RPA to second order in U. First note that if we sum the spin and 
charge susceptibilities over all wave vectors q and all Matsubara frequencies 
iujn,^ we obtain local, equal-time correlation functions, namely 

]^ XI XI X5p(q'^^") = ((% - ^if) = (^t) + i^l) - 2 ("T"i) (1-7) 

q iuJn 

and 

q iuJn 

(1.8) 

where on the right-hand side, we used the Pauli principle n^ = (cl^Ca) {clc„) = 
cJ.Ccr — c^clcaCcr = cl-Ccr = Ua that foUows from c^cl ~ CcrCa- — 0. This is the 
simplest version of the Pauli principle. Full antisymmetry is another mat- 
ter [31,32]. We call the first of the above displayed equations the local spin 
sum-rule and the second one the local charge sum-rule. For RPA, adding the 
two sum rules yields 

^ X X iXspifi,i^n) + Xch{q,ioJ„)) = (1.9) 

q icun 



T V- / xa{q) , xo{q) \ ^ 2 



N^\l-^Uxoiq) l + kUXoiq) 



2n-n\ (1.10) 



Since the non-interacting susceptibility Xo('?) satisfies the sum rule, we see by 
expanding the denominators that in the interacting case it is violated already 
to second order in U because Xo{q) being real and positive, (See Eq.(1.22)), 
the quantity ^„Xii{q)^ cannot vanish. 

How can we go about curing this violation of the Pauli principle while not 
damaging the conserving aspects? The simplest way is to proceed in the spirit 



^ In other references we often use ig„ instead of ia;„ to denote the Matsubara 
frequency corresponding to wave vector q. 
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of Fermi liquid theory and assume that the effective interaction (irreducible 
vertex in the jargon) is renormalized. This renormalization has to be different 
for spin and charge so that 

In practice x il) is the same^ as the Lindhard function xo{(l) Eq.(1.5) for 
U = but, strictly speaking, there is a constant self-energy term that is 
absorbed in the definition of /i [20]. We are almost done with the collective 
modes. Substituting the above expressions for XspIq) a-nd Xchil) in the two 
sum-rules, local-spin and local-charge appearing in Eqs.(1.7,1.8), we could 
determine both Usp and Uch if we knew (n-fn^) . The following ansatz 

U,p{n^){n^)=U{n^n^) (1.13) 

gives us the missing equation. Now notice that Usp, or equivalently (n-fni) 
depending on which of these variables you want to treat as independent, is 
determined self-consistently. That explains the name of the approach, "Two- 
Particle-Self-Consistent" . Since the the sum-rules are satisfied exactly, when 
we add them up the resulting equation, and hence the Pauli principle, will also 
be satisfied exactly. In other words, in Eq.(l.lO) that follows from the Pauli 
principle, we now have Usp and Uch on the left-hand side that arrange each 
other in such a way that there is no violation of the principle. In standard 
many-body theory, the Pauli principle (crossing symmetry) is achieved in a 
much more complicated way by solving parquet equations. [31,32] 

The ansatz Eq.(1.13) is inspired from the work of Singwi [33,34] and was 
also found independently by M. R. Hedeyati and G. Vignale [35]. The whole 
procedure is not as arbitrary as it may seem and we justify this in more detail 
in section 1.2.5 with the formal derivation. For now, let us just add a few 
physical considerations. 

Remark 1.1. Since Usp and Uch are renormalized with respect to the bare 
value, one might have expected that one should use the dressed Green's func- 
tions in the calculation of xo (q) ■ It is explained in appendix A of Ref. [7] 
that this would lead to a violation of the results X^(q = O^"^) — and 
X^(q = 0,w) = 0. In the present approach, the f-sum rule 

/ —^Xch.sp (q,^) = lim r^ (e— "" - e-"") iio„Xch,sp (q,*c.„) (1.14) 
= 1^ X! (*^k+q + ek-q - 2ek) rikCT (1-15) 

ko- 



^ The meaning of the superscripts differs from that in Ref. [7]. Superscripts (2) (1) 
here correspond respectively to (1) (0) in Ref. [7] 
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is satisfied with nj^^ — n^^, 

non-interacting case since it is computed from G^. ^ 



Remark 1.2. Usp{n^) {n^) — U (n^ni) can be understood as correcting tfie 
Hartree-Fock factorization to obtain the correct double occupancy. Expressing 
the irreducible vertex in terms of an equal-time correlation function is inspired 
by the approach of Singwi [33] to the electron gas. But TPSC is different since 
it also enforces the Pauli principle and connects to a local correlation function, 
namely (n^n^) . 

1.2.2 Mermin- Wagner, Kanamori-Brueckner and benchmarking 
spin and charge fluctuations 

The results that we found for spin and charge fluctuations have the RPA form 
but the renormalized interactions Usp and Uch rnust be computed from 

and 

With the ansatz Eq.(1.13), the above system of equations is closed and the 
Pauli principle is enforced. The first of the above equations is solved self- 
consistently with the Ugp ansatz. This gives the double occupancy (n-|-n^) 
that is then used to obtain Uch from the next equation. The fastest way to 
numerically compute x (?) is to use fast Fourier transforms [36]. 

These TPSC expressions for spin and charge fluctuations were obtained by 
enforcing the conservations laws and the Pauli principle. In particular, TPSC 
satisfles the f-sum rule Eq.(1.15). But we obtain for free a lot more, namely 
Kanamori-Brueckner renormalization and the Mermin- Wagner theorem. 

Let us begin with Kanamori-Brueckner renormalization of U . Many years 
ago, Kanamori in the context of the Hubbard model [2] , and Brueckner in the 
context of nuclear physics, introduced the notion that the bare U corresponds 
to computing the scattering of particles in the flrst Born approximation. In 
reality, we should use the full scattering cross section and the effective U 
should be smaller. From Kanamori's point of view, the two-body wave function 
can minimize the effect of U by becoming smaller to reduce the value of the 
probability that two electrons are on the same site. The maximum energy that 
this can cost is the bandwidth since that is the energy difference between a 



^ For the conductivity with vertex corrections [36] , the f-sum rule with riko- obtained 
from G"' is satisfied. 



Usr. - ^^r77 (1-19) 
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one-body wave function with no nodes and one with the niaxinium allowed 
number. Let us see how this physics comes out of our results. Far from phase 
transitions, we can expand the denominator of the local moment sum-rule 
equation to obtain 

|EE^'''('^)(l + ^^^P^^''(9))="-2^(n,)(n,). (1.18) 

Since ^ ^ '^iui^ X^^Kl) = n — 2 (n-f) (n^), we can solve for Usp and obtain * 

U 
1 + AU 

We see that at large U, Usp saturates to 1/yl, which in practice we find to 
be of the order of the bandwidth. For those that are familiar with diagrams, 
note that the Kanamori-Brueckner physics amounts to replacing each of the 
interactions U in the ladder or bubble sum for diagrams in the particle-hole 
channel by infinite ladder sums in the particle-particle channel [37]. This is 
not quite what we obtain here since (x^^') is in the particle-hole channel, 
but in the end, numerically, the results are close and the Physics seems to be 
the same. One cannot make strict comparisons between TPSC and diagrams 
since TPSC is non-perturbative. 

While Kanamori-Brueckner renormalization, or screening, is a quantum 
effect that occurs even far from phase transitions, when we are close we need 
to worry about the Mermin- Wagner theorem. To satisfy this theorem, approx- 
imate theories must prevent (S*^) from becoming infinite, which is equivalent 
to sloping {n^fi]^) from taking unphysical values. This quantity is positive 
and bounded by its value for U = oo and its value for non-interacting sys- 
tems, namely < {n^n{) < n^/4. Hence, the right-hand side of the local- 
moment sum-rule Eq.(1.16) is contained in the interval \n,n— -kn^] ■ To see 
how the Mermin- Wagner theorem is satisfied, write the self-consistency con- 
dition Eq.(1.16) in the form 

n-2{n^n^). (1.21) 



^V 


x'^'Hq) 


iV^l- 


Itt ("T"4> ,,1-11 („] 



Consider increasing {n^rn) on the left-hand side of this equation. The denom- 
inator becomes smaller, hence the integral larger. To become larger, (n^ni) 
has to decrease on the right-hand side. There is thus negative feedback in 
this equation that will make the self-consistent solution finite. This, however, 
does not prevent the expected phase transition in three dimensions [38]. To 

* There is a misprint of a factor of 2 in Ref. [7]. It is corrected in Rcf. [28]. 
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see this, wo need to look in more details at the phase space for the integral in 
the sum rule. 

As we know from the spectral representation for %, 

, . , fdu' Xch.sp (q^'^') f du' w'x"h,.p (q^^') ,, ^^, 

J TT UJ' -tUJn J TT (cj')^ + (w„)^ 

the zero Matsubara frequency contribution is always the largest. There, we 
find the so-called Ornstein-Zernicke form for the susceptibility. 

Ornstein-Zernicke form Let us expand the denominator near the point where 
1— ^Uspx'^^HQfi) — 0- The wave vector Q is that where x*-^-* is maximum. 



We find [38], 



x^'HQfi) 



2'~^spX' 4l^sp gQ2 </ 2^sp Q^ ^ 

e 



i^f2 2 ■ / ' (1-23) 

where all functions and derivatives in the denominator are evaluated at 
(Q,0) and where, on dimensional grounds. 



—^Usp -^^^7' 7 [I ^u,px'-'KQ,o)) (1-24) 

scales (noted ~) as the square of a length, ^, the correlation length. That 
length is determined self-consistently. Since LUgp ^ ^^^, all finite Mat- 
subara frequency contributions are negligible if 2ttT/u!sp ^ 2i:TS,^ ^ 1. 
That condition in the form ujgp ^ T justifies the name of the regime we 
are interested in, namely the renormalized classical regime. The classical 
regime of a harmonic oscillator occurs when uj <^ T. The regime here is 
"renormalized" classical because at temperatures above the degeneracy 
temperature, the system is a free classical gas. As temperature decreases 
below the Fermi energy, it becomes quantum mechanical, then close to 
the phase transition, it becomes classical again. 

Substituting the Ornstein-Zernicke form for the zero Matsubara frequency 
susceptibility in the self-consistency relation Eq.(1.16), we obtain 

d'^q 1 

(2nfq^+t 



T I .. ^, . , .-2 =g (1-25) 



where C contains all non-zero Matsubara frequency contributions as well as 
n — 2 {n^fi]^) . Since C is finite, this means that in two dimensions {d ~ 2), 
it is impossible to have ^^^ = on the left-hand side otherwise the integral 
would diverge logarithmically. This is clearly a dimension-dependent state- 
ment that proves the Mermin-Wagner theorem. In two-dimensions, we see 
that the integral gives a logarithm that leads to 
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Fig. 1.1. Wave vector (q) dependence of the spin and charge structure factors for 
different sets of parameters. Solid lines are from TPSC and symbols are QMC data. 
Monte Carlo data for n = 1 and U = St are for 6x6 clusters and T = 0.5i; all other 
data are for 8x8 clusters and T = 0.2t. Error bars are shown only when significant. 
From Ref. [6]. 



e ^ exp (C'/T) . 

where in general, C can be temperature dependent [38]. When C is not 
temperature dependent, the above result is similar to what is found at strong 
coupling in the non-linear sigma model. The above dimensional analysis is a bit 
expeditive. A more careful analysis [39,40] yields prefactors in the temperature 
dependence of the correlation length. 

The set of TPSC equations for spin and charge fluctuations Eqs. (1.16, 1.17, 1.13) 
is rather intuitive and simple. The agreement of calculations with benchmark 
QMC calculations is rather spectacular, as shown in Fig. (1.1). There, one can 
see the results of QMC calculations of the structure factors, i.e. the Fourier 
transform of the equal-time charge and spin correlation functions, compared 
with the corresponding TPSC results. This figure allows one to watch the 
Pauli principle in action. At t/ = At, Fig. (1.1a) shows that the charge structure 
factor does not have a monotonic dependence on density. This is because, as 
we approach half-filling, the spin fluctuations are becoming so large that the 
charge fluctuations have to decrease so that the sum still satisfies the Pauli 
principle, as expressed by Eq.(l.lO). 
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More comparisons may be found in Refs. [29] and [6,7,41,42]. This kind 
of agreement is found even at couplings of the order of the bandwidth and 
when second-neighbor hopping t' is present [43,44]. 

Remark 1.3. Even though the entry in the renormalized classical regime is 
well described by TPSC [45], equation (1.13) for Usp fails deep in that regime 
because Z"^^^ becomes too different from the true self-energy. At n = 1, t' = 
0, deep in the renormalized classical regime, Usp becomes arbitrarily small, 
which is clearly unphysical. However, by assuming that {n^n\^) is temperature 
independent below Tx^ a property that can be verified from QMC calculations, 
one obtains a qualitatively correct description of the renormalized-classical 
regime. One can even drop the ansatz and take (n-|-n^) from QMC on the 
right-hand side of the local moment sum-rule Eq.(1.16) to obtain Usp- 

1.2.3 Self-energy 

Collective charge and spin excitations can be obtained accurately from Green's 
functions that contain a simple self-energy, as we have just seen. Such modes 
are determined more by conservations laws than by details of the self-energy, 
especially at finite temperature where the lowest fermionic Matsubara fre- 
quency is not zero. The self-energy on the other hand is much more sensitive 
to collective modes since these are important at low frequency. The second 
step of TPSC is thus to find a better approximation for the self-energy. This 
is similar in spirit to what is done in the electron gas [8] where plasmons are 
found with non-interacting particles and then used to compute an improved 
approximation for the self-energy. This two step process is also analogous to 
renormalization group calculations where renormalized interactions are eval- 
uated to one-loop order and quasiparticle renormalization appears only to 
two-loop order [46-48]. 

The method to derive the result is justified using the formal derivation [20] 
that appears in Sect. 1.2. 5. If you are familiar with diagrams, you can under- 
stand physically the result by looking at Fig. 1.2 that shows the exact diagram- 
matic expressions for the three-point vertex (green triangle) and self-energy 
(blue circle) in terms of Green's functions (solid black lines) and irreducible 
vertices (red boxes). The bare interaction U is the dashed line. One should 
keep in mind that we are not using perturbation theory despite the fact that 
we draw diagrams. Even within an exact approach, the quantities defined in 
the figure have well defined meanings. The numbers on the figure refer to spin, 
space and imaginary time coordinates. When there is an over-bar, there is a 
sum over spin and spatial indices and an integral over imaginary time. 

In TPSC, the irreducible vertices in the first line of Fig. 1.2 are local, i.e. 
completely momentum and frequency independent. They are given by Usp and 
Uch- If we set point 3 to be the same as point 1, then we can obtain directly the 
TPSC spin and charge susceptibilities from that first line. In the second line of 
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Fig. 1.2. Exact expression for the three point vertex (green triangle) in the first line 
and for the self-energy in the second line. Irreducible vertices are the red boxes and 
Green's functions solid black lines. The numbers refer to spin, space and imaginary 
time coordinates. Symbols with an over-bard are summed/integrated over. The self- 
energy is the blue circle and the bare interaction U the dashed line. 
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Fig. 1.3. Exact self-energy in terms of the Hartree-Fock contribution and of the 
fully reducible vertex F represented by a textured box. 



the figure, the exact expression for the self-energy is displayed^. The first term 
on the right-hand side is the Hartree-Fock contribution. In the second term, 
one recognizes the bare interaction U at one vertex that excites a collective 
mode represented by the green triangle and the two Green's functions. The 
other vertex is dressed, as expected. In the electron gas, the collective mode 
would be the plasmon. If we replace the irreducible vertex using Ugp and Uch 
found for the collective modes, we find that here, both types of modes, spin 
and charge, contribute to the self-energy [41]. 

There is, however, an ambiguity in obtaining the self-energy formula [49]. 
Within the assumption that only Usp and Uch enter as irreducible particle-hole 
vertices, the self-energy expression in the transverse spin fiuctuation channel 

^ In the Hubbard model the Fock term cancels with the same-spin Hartree term 
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is different. What do we mean by that? Consider the exact formula for the seff- 
energy represented symboUcahy by the diagram of Fig. 1.3. In this figure, the 
textured box is the fuUy reducible vertex F {q, k — k' , k + k' — q) that depends 
in general on three momentum-frequency indices. The longitudinal version of 
the self-energy corresponds to expanding the fully reducible vertex in terms 
of diagrams that are irreducible in the longitudinal (parallel spins) channel 
illustrated in Fig. 1.2. This takes good care of the singularity of F when its first 
argument q is near (tt, tt) . The transverse version [20, 49] does the same for 
the dependence on the second argument k — k', which corresponds to the other 
(antiparallel spins) particle-hole channel. But the fully reducible vertex obeys 
crossing symmetry. In other words, interchanging two fermions just leads to 
a minus sign. One then expects that averaging the two possibilities gives a 
better approximation for F since it preserves crossing symmetry in the two 
particle- hole channels [49]. By considering both particle- hole channels only, 
we neglect the dependence oi F on k + k' — q because the particle-particle 
channel is not singular. The final formula that we obtain is [49] 

i:(2)(fc) = [/n_, + f ^ E i^UspXspiq) + UehXchiq)] G^^\k + q), (1.26) 

where n^^ is the average single-spin occupation. The superscript (2) reminds 
us that we are at the second level of approximation. Ga is the same Green's 
function as that used to compute the susceptibilities x^'^K'l)- Since the self- 
energy is constant at that first level of approximation, this means that Go- 
is the non-interacting Green's function with the chemical potential that gives 
the correct filling. That chemical potential fi^^'' is slightly different from the 
one that we must use in (G^^') ~ za;„ -I- /i*^^^ — Sk — -S'*-^-' to obtain the 
same density [50]. Estimates of /i*^^' may be found in Ref. [20,50]). Further 
justifications for the above formula are given below in Sect. 1.2. 4. 

But before we come up with more formalism, we check that the above 
formula is accurate by comparing in Fig. 1.4 the spectral weight (imaginary 
part of the Green's function) obtained from Eq.(1.26) with that obtained from 
Quantum Monte Carlo calculations. The latter are exact within statistical 
accuracy and can be considered as benchmarks. The meaning of the curves 
are detailed in the caption. The comparison is for half-filling in a regime 
where the simulations can be done at very low temperature and where a non- 
trivial phenomenon, the pseudogap, appears. This all important phenomenon 
is discussed further below in subsection 1.2.6 and in the first case study. Sect. 
1.3.1. In the third panel, we show the results of another popular Many-Body 
Approach, the FLuctuation Exchange Approximation (FLEX) [15]. It misses 
[51] the physics of the pseudogap in the single-particle spectral weight because 
it uses fully dressed Green's functions and assumes that Migdal's theorem 
applies, i.e. that the vertex does not need to be renormalized consequently 
Ref. [7, 52] . The same problem exists in the corresponding version of the GW 
approximation. [53] 
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Fig. 1.4. Single-particle spectral weight A(k, a;) for U = 4, jB — 5, n — 1, and all 
independent wave vectors Ic of an 8 x 8 lattice. Results obtained from maximum 
entropy inversion of Quantum Monte Carlo data on the left panel, from TPSC in 
the middle panel and form the FLEX approximation on the right panel. (Relative 
error in all cases is about 0.3%). Figure from Ref. [49] 



Remark 1.4- The dressing of one vertex in the second line of Fig. 1.2 means 
that we do not assume a Migdal theorem. Migdal's theorem arises in the case 
of electron-phonon interactions [54]. There, the small ratio m/M, where m 
is the electronic mass and M the ionic mass, allows one to show that the 
vertex corrections are negligible. This is extremely useful in formulating the 
Eliashberg theory of superconductivity. 

Remark 1.5. In Refs. [7,49] we used the notation E'^^^ instead of S'^'^\ The 
notation of the present paper is the same as that of Ref. [20] 



1.2.4 Internal accuracy checks 

How can we make sure that TPSC is accurate? We have shown sample compar- 
isons with benchmark Quantum Monte Carlo calculations, but we can check 
the accuracy in other ways. For example, we have already mentioned that the 
f-sum rule Eq.(1.15) is exactly satisfied at the first level of approximation (i.e. 
with nj^ on the right-hand side) . Suppose that on the right-hand side of that 
equation, one uses nt obtained from G^"^' instead of the Fermi function. One 
should find that the result does not change by more than a few percent. This 
is what happens when agreement with QMC is good. 

When we are in the Fermi liquid regime, another way to verify the accuracy 
of the approach is to verify if the Fermi surface obtained from G^"^^ satisfies 
Luttinger's theorem very closely. 
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Finally, there is a consistency relation between one- and two-particle quan- 
tities {S and (n-^nj^)). The relation 

§ E E ^(k, ^w„)G(k, zc.„)e-^""°" = ^Tr (SG) - U {n^n:,) (1.27) 

k n 

should be satisfied exactly for the Hubbard model. This result follows from 
the definition of self-energy and is derived in Eq.(1.40) below. In standard 
many-body books [54] , it is encountered in the calculation of the free energy 
through a coupling-constant integration. In TPSC, it is not difficult ^° to show 
that the following equation 

iTr(i:(2)G(i))=C/(n^n;) (1.28) 

is satisfied exactly with the self-consistent U {n^n^) obtained with the sus- 
ceptibilities^^. An internal accuracy check consists in verifying by how much 
^Tr (i^^^^G^^)) differs from ^Tr (I7(^)G(-^^) . Again, in regimes where we have 
agreement with Quantum Monte Carlo calculations, the difference is only a 
few percent. 

The above relation between S and (n-^n^) gives us another way to justify 
our expression for S^'^' . Suppose one starts from Fig. 1.2 to obtain a self- 
energy expression that contains only the longitudinal spin fluctuations and 
the charge fluctuations, as was done in the first papers on TPSC [6]. One 
finds that each of these separately contributes an amount U {n-^n^} /2 to the 
consistency relation Eq.(1.28). Similarly, if we work only in the transverse spin 
channel [20,49] we find that each of the two transverse spin components also 
contributes U (n^n^) /2 to ^Tr (Z'^^^G'-^') . Hence, averaging the two expres- 
sions also preserves rotational invariance. In addition, one verifies numerically 
that the exact sum rule (Ref. [7] Appendix A) 

J, / 

— S':" (k,c.') = C/2n_, (1 - n_.) (1.29) 

TT 

determining the high-frequency behavior is satisfied to a higher degree of 
accuracy with the symmetrized self-energy expression Eq. (1.26). 

Eq. (1.26) for S^^'> is different from so-called Berk-Schrieffer type expres- 
sions [56] that do not satisfy ^■^ the consistency condition between one- and 
two-particle properties, ^^r (SG) = U {n^n^) . 

Remark 1.6. Schemes, such as FLEX, that use on the right-hand side G^^^ 
are thermodynamically consistent (Sect. 1.4.4) and might look better. How- 
ever, as we just saw, in Fig. 1.4, FLEX misses some important physics. The 



^" Appendix B or Ref. [7] 

^^ FLEX does not satisfy this consistency requirement. See Appendix E of [7]. In 
fact double-occupancy obtained from EG can even become negative [55]. 
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(See Ref. [7] Appendix E) 
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reason [7] is that the vertex entering the self-energy in FLEX is not at the 
same level of approximation as the Green's functions. Indeed, since the latter 
contain self-energies that are strongly momentum and frequency dependent, 
the irreducible vertices that can be derived from these self-energies should 
also be frequency and momentum dependent, but they are not. In fact they 
are the bare vertices. It is as if the quasi-particles had a lifetime while at the 
same time interacting with each other with the bare interaction. Using dressed 
Green's functions in the susceptibilities with momentum and frequency inde- 
pendent vertices leads to problems as well. For example, the conservation law 
Xsp,ch (q = 0,iw„) = is violated in that case, as shown in Appendix A of 
Ref . [7] . Further criticism of conserving approaches appears in Appendix E of 
Ref. [7] and in Ref. [20]. 

1.2.5 A more formal derivation 

Details of a more formal derivation may be found in Ref. [57]. For complete- 
ness we repeat some of the derivation. The reader more interested in the 
physics may skip that section. The first two subsections present some general 
formalism that is then used in the following two subsections to derive TPSC. 

Single-particle properties 

Following functional methods of the Schwinger school [30,58,59], we begin with 
the generating function with source fields (pu and field destruction operators 
?/) in the grand canonical ensemble 



In Z \d)\ = In Tr 



3-^(^-^^)t, (e-'^-(i)Mi'2)M2)\l (^30) 



We adopt the convention that 1 stands for the position and imaginary time 
indices (ri,Ti). The over-bar means summation over every lattice site and 
integration over imaginary-time from to /3, and a summation over spins. T,- 
is the time-ordering operator. 

The propagator in the presence of the source field is obtained from func- 
tional differentiation 



G. (1, 2; {<t>}) = - (^. (1) i^l (2))^ = -^^§^- (1-31) 

From now on, the time- ordering operator in averages, (), is implicit. Physically 
relevant correlation functions are obtained for {</>} = but it is extremely 
convenient to keep finite {c/>} in intermediate steps of the calculation. 

Using the equation of motion for the field ?/; and the definition of the self- 
energy, one obtains the Dyson equation in the presence of the source field [60] 

(G(7^ - (/)) G = 1 + Z-G ; G-^ ^ Gq^ - 4) - S (1.32) 



18 A.-M.S. Tremblay 

where, from the commutator of the interacting part of the Hubbard Hamilto- 
nian _ff, one obtains 

S, (1, T; {4>}) G, (T, 2; {<(>}) - -U (^^, (l+) ^_, (1) ^„ (1) ^t (2)^^ . 

(1.33) 
The imaginary time in 1+ is infinitesimally larger than in 1. 

Response functions 

Response (four-point) functions for spin and charge excitations can be ob- 
tained from functional derivatives {SG/S(f)) of the source-dependent propaga- 
tor. Following the standard approach and using matrix notation to abbreviate 
the summations and integrations we have, 

GG-^ = 1 (1.34) 

Using the Dyson equation (1.32) G^^ — Gq ^ — <p — S this may be rewritten 



^-^ = -G^-^G = G^G + G^-^G, (1.36) 

ocp ocp oq) 

where the symbol - in G-G reminds us that the neighboring labels of the 
propagators have to be the same as those of the <j) in the functional derivative. 
If perturbation theory converges, we may write the self-energy as a functional 
of the propagator. From the chain rule, one then obtains an integral equation 
for the response function in the particle-hole channel that is the analog of the 
Bethe-Salpeter equation in the particle-particle channel 



-. , — Gr" Gr -|- O 

0(p 



SSSG 

JgJ^) 



G. (1.37) 



The labels of the propagators in the last term are attached to the self energy, 
as in Eq.(1.36) ^^ . Vertices appropriate for spin and charge responses are given, 
respectively, by 






^'P^6G, 6G, ' ^^''"<5G, +(5G/ ^^'^^^ 



To remind ourselves of this, we may also adopt an additional vertical matrix 



notation convention and write Eq.(7) as j^ = G~G + G 



5_S 
SG 
~STT 

5<p 



G. 
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TPSC First step: two-particle self-consistency for G'^-'^^ X'(^\ 

rw ^ Usp and r^l^ = [7,^ 

In conserving approximations, the self-energy is obtained from a functional 
derivative S [G] = S<P [G] /SG of ^ the Luttinger-Ward functional, which is 
itself computed from a set of diagrams. To liberate ourselves from diagrams, 
we start instead from the exact expression for the self-energy, Eq.(1.33) and 
notice that when label 2 equals l"*", the right-hand side of this equation is 
equal to double-occupancy (n^ni). Factoring as in Hartree-Fock amounts to 
assuming no correlations. Instead, we should insist that {n^n^) be obtained 
self-consistently. After all, in the Hubbard model, there are only two local 
four point functions: {n^n\) and {n^n^) = (n^n^) . The latter is given exactly, 
through the Pauli principle, by (n^n^) = {n^n^) = (n^) = (n^) = n/2, when 
the filling n is known. In a way, (n-fni) in the self-energy equation (1.33), 
can be considered as an initial condition for the four point function when 
one of the points, 2, separates from all the others which are at 1. When that 
label 2 does not coincide with 1, it becomes more reasonable to factor a la 
Hartree-Fock. These physical ideas are implemented by postulating 



41) (1, 1; m) G(i) (1, 2; {0}) = A^^jG^ (l, 1+; {</-}) G^^) (1, 2; {0}) 

(1.39a) 

where ^{0} depends on external field and is chosen such that the exact result 

14 

S, (1, T; {^}) G„ (T, 1+; {(/)}) = U {n^ (1) n^ (1))^ (1.40) 

is satisfied. It is easy to see that the solution is 

^''' - ''(%(i)),(n,(i))/ ^'-''^ 

Substituting ^{0} back into our ansatz Eq.(1.13) wc obtain our first approx- 
imation for the self-energy by right-multiplying by I Gij 



4') (1, 2; {0}) = A{^}GL'1 (1, 1+; {0}) -5 (1 - 2) . (1.42) 

We are now ready to obtain irreducible vertices using the prescription 
of the previous section, Eq.(1.38), namely through functional derivatives of 
E with respect to G. In the calculation of Ugp, the functional derivative of 
(ni-ni) I ((%) (n^)) drops out, so we are left with ^^ 



^"^ See footnote (14) of Ref. [20] for a discussion of the choice of limit 1+ vs 1 . 
^^ For n > 1, all particle occupation numbers must be replaced by hole occupation 



numbers. 
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U - —^ 



5S['^ 



{4>}=o t 



= A{,j.o = C/7^^. (1.43) 



The renormalization of this irreducible vertex may be physically understood as 
coming from Kanamori-Brucckncr screening [7] . This completes the derivation 
of the ansatz that was missing in our first derivation in section 1.2.1. 

The functional-derivative procedure generates an expression for the charge 
vertex Uch which involves the functional derivative of {n-^ni) / {{n^} (%)) 
which contains six point functions that one does not really know how to 
evaluate. But, if we again assume that the vertex Uch is a constant, it is sim- 
ply determined by the requirement that charge fluctuations also satisfy the 
fluctuation-dissipation theorem and the Pauli principle, as in Eq.(1.17). 

Note that, in principle, S^^' also depends on double-occupancy, but since 
S^^' is a constant, it is absorbed in the definition of the chemical potential and 
we do not need to worry about it in this case. That is why the non-interacting 
irreducible susceptibility x (?) = Xo{q) appears in the expressions for the 
susceptibility, even though it should be evaluated with G^^-* that contains E^^\ 
A rough estimate of the renormalized chemical potential (or equivalently of 
S^^'^), is given in the appendix of Ref. [20] . One can check that spin and charge 
conservation are satisfied by our susceptibilities. 

TPSC Second step: an improved self-energy S^^^ 

Collective modes are emergent objects that are less influenced by details of 
the single-particle properties than the other way around. We thus wish now 
to obtain an improved approximation for the self-energy that takes advantage 
of the fact that we have found accurate approximations for the low-frequency 
spin and charge fluctuations. We begin from the general definition of the self- 
energy Eq.(1.33) obtained from Dyson's equation. The right-hand side of that 
equation can be obtained either from a functional derivative with respect to an 
external field that is diagonal in spin, as in our generating function Eq.(1.30), 
or by a functional derivative of (V'-o- (1) V'J (2)) j, with respect to a transverse 
external field (j)t- 

Working first in the longitudinal channel, the right-hand side of the general 
definition of the self-energy Eq.(1.33) may be written as 



r, (1,1) G, (1,2) =-(7 



<5G, (1,2; {,/.}) 



(1+,1) 



-G_4l,l+)G,(1,2) 

(1.44) 



The last term is the Hartree-Fock contribution. It gives the exact result for the 
self-energy in the limit lu — > oo. [7] The 5Ga / 5(f)^„ term is thus a contribution 
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to lower frequencies and it conies from the spin and charge fluctuations. Right- 
multiplying the last equation by G"^ and replacing the lower energy part 
5Ga/54>-(j by its general expression in terms of irreducible vertices, Eq.(1.37) 
we find 

i:(2)(l,2) = [/GL^^(f,l+)<5(l-2) (1.45) 



-t/GW(l.3) 



(5^(^'(3,2;{0}) 



5G^^^(4,5;W) 



X 



<5G«(4,5;{</-}) 



(.HO ^^-^(^"'1) 



{0}=O_ 



Every quantity appearing on the right-hand side of that equation has been 
taken from the TPSC results. This means in particular that the irreducible 
vertices SSa /SG)^, are at the same level of approximation as the Green func- 
tions Gct and self-energies Sa ■ In approaches that assume that Migdal's 
theorem applies to spin and charge fluctuations, one often sees renormalized 
Green functions G*-^' appearing on the right-hand side along with unrenor- 
malizcd vertices, SSa/SGai — > U. In terms of Uap and Uch in Fourier space, 
the above formula [41] reads. 



(1.46) 

The approach to obtain a self-energy formula that takes into account both 
longitudinal and transverse fluctuations is detailed in Ref . [20] . Crossing sym- 
metry, rotational symmetry and sum rules and comparisons with QMC dictate 
the final formula for the improved self-energy S^'^' as we have explained in 
Sect. (1.2. 3). 

1.2.6 Pseudogap in the renormalized classical regime 

When we compared TPSC with Quantum Monte Carlo simulations and with 
FLEX in Fig. 1.4 above, perhaps you noticed that at the Fermi surface, the 
frequency dependent spectral weight has two peaks instead of one. In addition, 
at zero frequency, it has a minimum instead of a maximum. That is called a 
pseudogap. A cartoon explanation [29] of this pseudogap is given in Fig. 1.5. At 
high temperature we start from a Fermi liquid, as illustrated in panel I. Now, 
suppose the ground state has long-range antiferromagnetic order as in panel 
III, in other words at a filling between half- filling and ric. In the mean- field 
approximation we have a gap and the Bogoliubov transformation from fermion 
creation-annihilation operators to quasi-particles has weight at both positive 
and negative energies. In two dimensions, because of the Mermin- Wagner 
theorem, as soon as we raise the temperature above zero, long-range order 
disappears, but the antiferromagnetic correlation length ^ remains large so we 
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Fig. 1.5. Cartoon explanation of the pseudogap due to precursors of long-range 
order. When the antiferromagnetic correlation length ^ becomes larger than the 
thermal de Broglie wavelength, there appears precursors of the T = Bogoliubov 
quasi-particles for the long-range ordered antiferromagnet. This can occur only in 
the renormalized classical regime, below the dashed line on the left of the figure. 



obtain the pseudogap illustrated in panel II. As we will explain analytically 
below, the pseudogap survives as long as ^ is much larger than the thermal 
de Broglie wave length S^th = vf/{t^T) in our usual units. At the crossover 
temperature Tx, the relative size of ^ and ^th changes and we recover the 
Fermi liquid. 

We now proceed to sketch analytically where these results come from 
starting from finite T. Details and more complete formulae may be found 
in Refs. [6,7,41,61]^^ We begin from the TPSC expression (1.26) for the self- 
energy. Normally one has to do the sum over bosonic Matsubara frequencies 
first, but the zero Matsubara frequency contribution has the correct asymp- 
totic behavior in fermionic frequencies «a;„ so that, as in Sect. 1.2. 2, one can 
once more isolate on the right-hand side the contribution from the zero Mat- 
subara frequency. In the renormalized classical regime then, we have 



U{kF,lCJn)^T f q"^-^ 



dq- 



1 



1 



q- +(, '^iuJn- Ekf+Q+q 



(1.47) 



where Q is the wave vector of the instability. ^^ This integral can be done 
analytically in two dimensions [7, 64] . But it is more useful to analyze limiting 
cases [41]. Expanding around the points known as hot spots where ek^+Q — Oj 
we find after analytical continuation that the imaginary part of the retarded 
self-energy at zero frequency takes the form 



^^ Note also the following study from zero temperature [62] 
^'^ This formula is similar to one that appeared in Ref. [63] 
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i:"«(kj.,0) ex -ttT [ d^-\^dq\\ , \^.2 ^{v'f%) (1-48) 

oc ^f -'^. (T49) 

Vp 

In the last line, we just used dimensional analysis to do the integral. 

The importance of dimension comes out clearly [41]. In d = 4, 17"^ (k^?, 0) 
vanishes as temperature decreases, d = 3 is the marginal dimension and in 
d = 2 we have that i7"^(ki?, 0) ex £,/£,th that diverges at zero temperature. In 
a Fermi liquid the quantity Z'"^(ki?,0) vanishes at zero temperature, hence 
in three or four dimensions one recovers the Fermi liquid (or close to one 
m d — 3). But in two dimensions, a diverging Z'"^(ki?,0) corresponds to a 
vanishingly small A{\<ip^oj = 0) as we can see from 

^('^'") = 7 r'n^'^'^^^^^.n,^ ,2 - (1-50) 

Fig. 31 of Ref. [29] illustrates graphically the relationship between the location 
of the pseudogap and large scattering rates at the Fermi surface. At stronger 
U the scattering rate is large over a broader region, leading to a depletion of 
ACkjOj) over a broader range of k values. 

Remark 1.7. Note that the condition £,/£,th ^ 1, necessary to obtain a large 
scattering rate, is in general harder to satisfy than the condition that cor- 
responds to being in the renormalized classical regime. Indeed, £,/S,th ^ 1 
corresponds T/vp 3> ^~^ while the condition LOgp ^ T for the renormal- 
ized classical regime corresponds to T ^ C^^, with appropriate scale factors, 
because Wgp scales as ^^^ as we saw in Eq. (1.23) and below. 

To understand the splitting into two peaks seen in Figs. 1.4 and 1.5 con- 
sider the singular renormalized contribution coming from the spin fluctuations 
in Eq. (1.47) at frequencies uj ^ vp^^^. Taking into account that contribu- 
tions to the integral come mostly from a region q < ^^^, one finds 



S"'{kp,co)=(Tjq^-'dq ^ ^ ^ 



(j2 +^ ^ J ik-a -eki=.+Q 

Z\2 



^ — £ki=.+Q 



(1.51) 



which, when substituted in the expression for the spectral weight (1.50) leads 
to large contributions when 

w-£k =0 (1.52) 

^ - Ekp+Q 

or, equivalently, 
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_ (£k + Ek^+q) ± V(£k -Ek^+q)^ +4Z\2 



(1.53) 



which, at (x) = 0, corresponds to the position of the hot spots^^. At finite 
frequencies, this turns into the dispersion relation for the antiferroniagnet [66]. 

It is important to understand that analogous arguments hold for any fluc- 
tuation that becomes soft because of the Mermin- Wagner theorem, [7, 67] in- 
cluding superconducting ones [7,50,57]. The wave vector Q would be different 
in each case. 

To understand better when Fermi liquid theory is valid and when it is 
replaced by the pseudogap instead, it is useful to perform the calculations 
that lead to Z"'^(ki?,0) ex ^/^th in the real frequency formalism. The details 
may be found in Appendix D of Ref. [7] . 



1.3 Case studies 

In this short pedagogical review it is impossible to cover all topics in depth. 
This section will nevertheless expand a bit on two important contributions 
of TPSC to problems of current interest, namely the pseudogap of cuprate 
superconductors and superconductivity induced by antiferromagnetic fluctu- 
ations. 

1.3.1 Pseudogap in electron-doped cuprates 

High-temperature superconductors are made of layers of Cu02 planes. The 
rest of the structure is commonly considered as providing either electron or 
hole doping of these planes depending on chemistry. At half-filling, or zero- 
doping, the ground state is an antiferromagnet. As one dopes the planes, one 
reaches a doping, so-called optimal doping, where the superconducting tran- 
sition temperature T^ is maximum. Let us start from optimal hole or elec- 
tron doping and decrease doping towards half-filling. That is the underdoped 
regime. In that regime, one observes a curious phenomenon, the pseudogap. 
What this means is that as temperature decreases, physical quantities behave 
as if the density of states near the Fermi level were decreasing. Finding an 
explanation for this phenomenon has been one of the major challenges of the 
field [68,69]. 

To make progress, we need a microscopic model for high-temperature su- 
perconductors. Band structure calculations [70, 71] reveal that a single band 
crosses the Fermi level. Hence, it is a common assumption that these mate- 
rials can be modeled by the one-band Hubbard model. Whether this is an 
oversimplification is still a subject of controversy [72-77]. Indeed, spectro- 
scopic studies [72, 78] show that hole doping occurs on the oxygen atoms. 



For comparisons with paramagnon theory see [65]. 
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The resulting hole behaves as a copper excitation because of Zhang-Rice [79] 
singlet formation. In addition, the phase diagram [80-85] and many proper- 
ties of the hole-doped cuprates can be described by the one-band Hubbard 
model. Typically, the band parameters that are used are: nearest-neighbor 
hopping t = 350 to 400 meV and next-nearest-neighbor hopping t' — —0.15 
to — 0.3i depending on the compound [70,71]. A third-nearest-neighbor hop- 
ping t" — — 0.5i' is sometimes added to fit finer details of the band struc- 
ture [71]. The second- neighbor hopping breaks particle-hole symmetry at the 
band structure level. 

In electron-doped cuprates, the doping occurs on the copper, hence there 
is little doubt that the single-band Hubbard model is even a better starting 
point in this case. Band parameters [86] are similar to those of hole-doped 
cuprates. It is sometimes claimed that there is a pseudogap only in the hole- 
doped cuprates. The origin of the pseudogap is indeed probably different in the 
hole-doped cuprates. But even though the standard signature of a pseudogap 
is absent in nuclear magnetic resonance [87] (NMR) there is definitely a pseu- 
dogap in the electron-doped case as well [88] , as can be seen in optical conduc- 
tivity [89] and in Angle Resolved Photoemission Spectroscopy (ARPES) [90]. 
As we show in the rest of this section, in electron-doped cuprates strong evi- 
dence for the origin of the pseudogap is provided by detailed comparisons of 
TPSC with ARPES as well as by verification with neutron scattering [91] that 
the TPSC condition for a pseudogap, namely ^ > £_th, is satisfied. The latter 
length makes sense from weak to intermediate coupling when quasi-particles 
exist above the pseudogap temperature. In strong coupling, i.e. for values of U 
larger than that necessary for the Mott transition, there is evidence that there 
is another mechanism for the formation of a pseudogap. This is discussed at 
length in Refs. [92,93] ^^. The recent discovery [94] that at sufficiently large U 
there is a first order transition in the paramagnetic state between two kinds 
of metals, one of which is highly anomalous, gives a sharper meaning to what 
is meant by strong-coupling pseudogap. 

Let us come back to modeling of electron-doped cuprates. Evidence that 
these are less strongly coupled than their hole-doped counterparts comes from 
the fact that a) The value of the optical gap at half- filling, ~ 1.5 eV, is smaller 
than for hole doping, ^ 2.0 cV [95]. b) In a simple Thomas- Fermi picture, the 
screened interaction scales like dfx/dn. Quantum cluster calculations [92] show 
that dix/dn is smaller on the electron-doped side, hence U should be smaller. 

c) Mechanisms based on the exchange of antiferromagnetic fluctuations with 
U/t at weak to intermediate coupling [16,42] predict that the superconducting 
Tc increases with U/t. Hence Tc should decrease with increasing pressure 
in the simplest model where pressure increases hopping t while leaving U 
essentially unchanged. The opposite behavior, expected at strong coupling 
where J — At^ /U is relevant [85,96], is observed in the hole-doped cuprates. 

d) Finally and most importantly, we have shown detailed agreement between 

^® See also conclusion of Ref. [29] . 
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Fig. 1.6. On the left, results of TPSC calculations [29,97] at optimal doping, x = 
0.15, corresponding to filling 1.15, for t = 350 meV, t' = -0.175t, t" = 0.05i, 
U = 5.75t, T = 1/20. The left-most panel is the magnitude of the spectral weight 
times a Fermi function, A (k, lj) f (uj) at uj — 0, so-called momentum-distribution 
curve (MDC). Red (dark black) indicates larger value and purple (light grey) smaller 
value. The next panel is A (k, lj) / (uj) for a set of fixed k values along the Fermi 
surface. These are so-called energy-dispersion curves (EDC). The two panels to the 
right are the corresponding experimental results [90] for Nd2-i:Cea;Cu04. Dotted 
arrows show the correspondence between TPSC and experiment. 



TPSC calculations [29,93,97] and measurements such as ARPES [90,98], 
optical conductivity [89] and neutron [91] scattering. 

To illustrate the last point, consider Fig. 1.6 that compares TPSC calcu- 
lations with experimental results for ARPES. Apart from a tail in the exper- 
imental results, the agreement is striking. ^'^. In particular, if there were no 
interaction, the Fermi surface would be a line (red) on the momentum dis- 
tribution curve (MDC). Instead, it seems to disappear at symmetrical points 
displaced from {n/2, n/2) . These points, so-called hot spots, are linked by the 
wave vector (tt, tt) to other points on the Fermi surface. This is where the an- 
tiferromagnetic gap would open first if there were long-range order. The pull 
back of the weight from w = at the hot spots is close to the experimental 
value: 100 meV for the 15% doping shown, and 300 meV for 10% doping (not 
shown). More detailed ARPES spectra and comparisons with experiment are 
shown in Ref. [29] . The value of the temperature T* at which the pseudogap 
appears [97] is also close to that observed in optical spectroscopy [89]. In ad- 
dition, the size of the pseudogap is about ten times T* in the calculation as 
well as in the experiments. For optical spectroscopy, vertex corrections (see 
Sect. 1.4.5) have to be added to be more quantitative. Experimentally, the 
value of T* is about twice the antifcrromagnetic transition temperature up 



^° Such tails tend to disappear in more recent laser ARPES measurements on hole- 
doped compounds [99]. 
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to a; = 0.13. That can be obtained [97] by taking tz — 0.03i for hopping in 
the third direction. Recall that in strictly two dimensions, there is no long- 
range order. Antiferromagnetism appears on a much larger range of dopings 
for electron-doped than for hole-doped cuprates. 

These TPSC calculations have predicted the value of the pseudogap tem- 
perature at X = 0.13 before it was observed experimentally [98] by a group 
unaware of the theoretical prediction. In addition, the prediction that ^ should 
scale like ^th at the pseudogap temperature has been verified in neutron scat- 
tering experiments [91] in the range x — 0.04 to x = 0.15. At that doping, 
which corresponds to optimal doping, T* becomes of the order of 100 K, more 
than four times lower than at a; = 0.04. The antiferromagnetic correlation 
length ^ beyond optimal doping begins to decrease and violate the scaling 
of S, with ^(/j. In that doping range, T* and the superconducting transition 
temperature are close. Hence it is likely that there is interference between the 
two phenomena [100], an effect that has not yet been taken into account in 
TPSC. 

An important prediction that one should verify is that inelastic neutron 
scattering will find over-damped spin fluctuations in the pseudogap regime 
and that the characteristic spin fluctuation energy will be smaller than ksT 
whenever a pseudogap is present. Equality should occur above T* . 

Finally, note that the agreement found in Fig. 1.6 between ARPES and 
TPSC is for U ^ 6i. At smaller values of U the antiferromagnetic correlations 
are not strong enough to produce a pseudogap in that temperature range. For 
larger U, the weight near (7r/2, 7r/2) disappears, in disagreement with experi- 
ments. The same value of U is found for the same reasons in strong coupling 
calculations with Cluster Perturbation Theory (CPT) [92] and with slave bo- 
son methods [101]. Recent first principle calculations [102] find essentially the 
same value of U. In that approach, the value of U is fixed, whereas in TPSC it 
was necessary to increase U by about 10% moving towards half-filling to get 
the best agreement with experiment. In any case, it is quite satisfying that 
weak and strong coupling methods agree on the value of U for electron-doped 
cuprates. This value of U is very near the critical value for the Mott transi- 
tion at half- filling [103]. Hence, antiferromagnetic fluctuations at finite doping 
can be very well described by Slater-like physics (nesting) in electron-doped 
cuprates. 

For recent calculations including the effect of the third dimension on the 
pseudogap see [104]. Finally, note that the analog of the above mechanism 
for the pseudogap has also been seen in two-dimensional charge-density wave 
dichalcogenides [105]. 

1.3.2 d-wave superconductivity 

In the BCS theory of superconductivity, pairs of electrons form because of 
an effective attraction mediated by phonons. The pairs then condense in 
a coherent state. The suggestion that superconductivity could arise from 
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purely repulsive forces goes back to Kohn and Luttinger who showed that 
for pairs of sufficiently high angular momentum, the screened Coulomb in- 
teraction in an electron gas could be attractive [106]. Just before the dis- 
covery of high-temperature superconductors, an extension of that idea was 
proposed [107-109]. The suggestion was that antiferromagnetic fluctuations 
present in the Hubbard model could replace the phonons in BCS theory and 
lead to d— wave superconductivity. This is difficult to prove beyond any doubt 
since superconductivity in this case does not arise at the mean-field level. 
Mean-field on the Hubbard model gives antiferromagnetism near half-filling 
but not superconductivity. In the high-temperature superconductors, the sit- 
uation is made even more difficult because of Mott Physics. Nevertheless, 
the question is well posed and, as we just saw, Mott physics might be less 
important in the electron-doped superconductors. 

To investigate how the pairing susceptibility is influenced by antiferromag- 
netic fluctuations in TPSC we proceed as follows [42,110]. The reader who did 
not go through the formal section 1.2.5 may skip the next paragraph without 
loss of continuity to read the physical results below. The few equations that 
appear below give details that are missing in the literature. 

Some details of the derivation: We work in Nambu space and add an off- 
diagonal source fleld 6 and 6* in the generating function Eq.(1.30). The 
transverse spin fluctuations are included by working with four by four ma- 
trices. The pair susceptibility in the normal state can be obtained from 
the second functional derivative of the generating function with respect 
to the off-diagonal source field, evaluated at zero source field. [42] In more 
detail, one proceeds as in the formal derivation in Sect. 1.2.5. The ex- 
pressions for the spin and charge susceptibilities are not modified. Once 
the two-particle quantities have been found as above, the second step 
of the approach [7, 49] consists in improving the approximation for the 
single-particle self-energy by starting from the exact expression where the 
high-frequency Hartrce-Fock behavior is explicitly subtracted 

^(^'^)-<;o)f(^«-M3,2) (1.54) 

/0 0\ ^G(l,3) 

The bold face objects are matrices in Nambu space. To be able to express 
the right-hand side of the above equation in terms of irreducible vertices, 
susceptibilities and powers of G, one differentiates GG^^ = I to obtain 
{5G/5cj)) G-i = -G {5G-'^/5cj)) . With the help of Dyson's equation on 
the right-hand side of the last equation as well as the chain rule, one finds 
an expression where one can replace every term by their value at the first 
step, namely Ugp and Uch for the irreducible low-frequency vertices as 
well as G)j (fc + q) and Xsp{q) , Xch{q) ■ For the diagonal piece of the self- 
energy at the second step, one then obtains Eq.(1.26) above or equivalently 
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Eq.(3) of Rcf. [20,49] by considering both longitudinal and transverse 
channels and requiring crossing symmetry of the fully-reducible vertex in 
the two particle-hole channels as well as consistency with the sum-rule 
Tr(i;(2)G'(i)) = 2U {n-^n^) [7]. The off-diagonal piece of the exact self- 
energy Eq.(1.54) on the other hand, reads 

^f,'^-t/GW(4,3)'^5i|M (1.55) 



-UG^l^ (4, 4+) ,5 (4 - 5) - UG''1^ (4, 3' 



^^'542 (3,5 



%(4+,4) 

The pairing susceptibility mediated by spin fluctuations may now be com- 
puted from the derivative with respect to the source field 

Gl^i'(l,2)G(f (4,3) (1.56) 
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^22 (5,3) 



8=0 



with the irreducible vertex JZ'Jj /'5G]^2 obtained by functional differenti- 
ation of Eq.(1.55). Neglecting 5/5(1)^ {6S12 I^Gi2 ) , which represents the 

influence of spin fluctuations on the local piece of the irreducible vertex, 
and including the transverse component, we find for the d-wave pair sus- 
ceptibility, Xd, the expression that appears in Eq.(l) of Ref. [42]. The 

TPSC expression for the pair susceptibility Xd contains the bubble part 

(2) 
and the first term of what would be an infinite series if Z'Jg in the ir- 

f2) 

reducible vertex could be differentiated with respect to G\2 instead of 

Why should we trust the results for the d-wave susceptibility obtained for 
TPSC? Let us look again at benchmarks. Fig. 1.7(a) displays the d-wave sus- 
ceptibility obtained from QMC calculations shown as symbols and from TPSC 
as lines. Because of the sign problem, it is not practical to do the QMC calcu- 
lations at lower temperatures. Nevertheless, the temperatures are low enough 
that we see a non-trivial effect, the appearance of a maximum in susceptibility 
at finite doping and a substantial increase with decreasing temperatures. The 
agreement between QMC and TPSC is to within a few percent and improves 
for lower values of U . When the interaction strength reaches the intermediate 
coupling regime, U — 6, deviations of the order of 20% to 30% may occur 
but the qualitative dependence on temperature and doping remains accurate. 
In TPSC the pseudogap is the key ingredient that leads to a decrease in Xd 
in the underdoped regime. This is easy to understand since the pseudogap 
leaves fewer states for pairing at the Fermi level. Another way to say this is 
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Fig. 1.7. (Color online) From Ref. [42]. (a) Comparisons between the d^2_y2 sus- 
ceptibility obtained from QMC simulations and from the approach described in the 
present work. QMC error bars are smaller than the symbols. Analytical results are 
joined by solid lines. Both calculations are for [/ = 4, a 6 x 6 lattice, and four different 
temperatures. The case U — 0, P = 4 is shown for reference. The size dependence 
of the results is small at these temperatures. The inset compares QMC and FLEX 
at [/ = 4, /3 = 4. (b) Contributions from the bubble (DOS) represented by squares 
and vertex represented by circles, (c) Estimate of Tc using the Thouless criterion for 
[/ = 4 and (7 = 6, f' = t" = 0. 



that the strong inelastic scattering that leads to the pseudogap is pair break- 
ing. The inset shows that previous spin-fluctuation calculations (FLEX) in 
two dimensions [16, 111] deviate both qualitatively and quantitatively from 
the QMC results. More specifically, in the FLEX approach Xd does not show 
a pronounced maximum at finite doping. This is because, as we have shown in 
Fig. 1.4, in FLEX there is no pseudogap in the single-particle spectral weight 
at the Fermi surface [49,51]. 

In two-dimensions, superconductivity is in the Kosterlitz-Thouless univer- 
sality class. Vortex physics that is absent in TPSC is important to understand 
the precise value of the transition temperature. Nevertheless, a necessary con- 
dition for this transition to occur is that there is a higher temperature at 
which pairs form, a sort of mean-field Tc. That Tc can be obtained from the 
so-called Thouless criterion, i.e. from the temperature at which the d-wave 
susceptibility diverges. This divergence occurs because of growing vertex cor- 
rections. Since TPSC contains only the first term in what would be an infinite 
ladder sum, we take Tc as the temperature at which the bubble (that we call 
DOS) and first term of the series (that we call vertex) become equal. This is 
illustrated in Fig. 1.7(b). At (3 = 1/64, the doping range between the inter- 
section of the two curves is below Tc. The resulting Tc versus doping is shown 



in Fig. 1.7(c). In that calculation. 



«T"i/ 



is fixed at its value at T* . 



There is clearly an additional approximation in finding Tc with TPSC that 
goes beyond what can be checked with QMC calculations. How can we be sure 
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that this is correct? First there is consistency with other weak-coupling ap- 
proaches. For example, Ref. [112] has shown consistency with cutoff renormal- 
ization group technique [24] for competing ferromagnetic, antiferromagnetic 
and d-wave superconductivity. Second, there is consistency as well with Quan- 
tum Cluster Approaches that are best at strong-coupling. Indeed, an extensive 
study as a function of system size by the group of Jarrell [81] has shown that 
for 10% doping, Tc is in the range 0.02, not far from 0.03 that can be read 
from Fig. 1.7(c). 

One of the major theoretical questions in the field of high-temperature 
superconductivity has been, "Is there d-wave superconductivity in the two- 
dimensional Hubbard model"? TPSC has contributed to answer this ques- 
tion. One of the most discouraging early results was that QMC simulations 
showed that the pairing susceptibility was smaller at finite U than at U — 0. 
This is clearly seen in Fig. 1.7(a). TPSC allows us to understand why. At 
/3 = 4, the bubble largely dominates and the effect can only be pair break- 
ing because of the inelastic scattering. The vertex, representing exchange of 
antiferromagnetic fluctuations analogous to the phonons in ordinary BCS the- 
ory, contributes 22% at most at zero doping and much less at larger doping. 
Nevertheless, it clearly increases the pair susceptibility to bring it in closer 
agreement with QMC. TPSC allows us to do the calculation at temperatures 
much lower than QMC and to verify that indeed the vertex eventually grows 
large enough to lead to a transition. We understand also that in this parameter 
range the dome shape comes from the fact that antiferromagnetic fluctuations 
can both increase pairing through the vertex and be detrimental through the 
pseudogap produced by the large self-energy. Antiferromagnetic fluctuations 
can both help and hinder d-wave superconductivity. 

Another question is whether the presence of a quantum critical point below 
the maximum of the superconducting dome plays a role in superconductivity. 
In the case we discussed above, long-range antiferromagnetic order appears 
at T = (not at finite T because of Mermin- Wagner) up to doping 6 = 0.17 
for L'^ = 4 [42] and 6 = 0.205 for t/ = 6 [36]. In this case, then, according to 
Fig. 1.7(c) the quantum critical point is far to the right of the maximum Tc 
but superconductivity can exist to the right of that point, the more so when 
U is larger. 

How general are the above results? This and many more questions on 
the conditions for magnetically mediated superconductivity were studied with 
TPSC on the square lattice at half-filling for second-neighbor hopping t' differ- 
ent from zero [110]. At t' = at half-filling, there is a pseudogap on the whole 
Fermi surface because of perfect nesting, so Tc vanishes. When t' is increased 
from zero, the pseudogap is not complete at half-filling and T^ is different from 
zero. In addition, for t' larger than 0.71, the Fermi surface topology changes 
and the dominant magnetic fluctuations are near (0, ±7r) , (±7r, 0) . 

Additional conclusions of the TPSC study of Ref. [110] as a function of 
t' and U at half-filling are as follows. First some qualitative conclusions that 
could be found from just the BCS gap equation with an interaction potential 
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given by the static component of the spin susceptibiUty [113]: the symmetry of 
the d-wave order parameter is determined by the wave vector of the magnetic 
fluctuations. Those that are near {tt,tt) lead to dx2_y2-wa,ve (Big) supercon- 
ductivity while those that are near (0,7r) induce d^y-wave (-B23) superconduc- 
tivity. The dominant wave vector for magnetic fluctuations is determined by 
the shape of the Fermi surface so dx2_y2-wave superconductivity occurs for 
values of t' that are relatively small while dxy-'wstve superconductivity occurs 
for i' > 1. Second, the maximum value that T^ can take as a function of t' 
increases with interaction strength. TPSC cannot reach the strong-coupling 
regime where Tc should decrease with U. 

One also finds [110] that, contrary to what is expected from BCS, the non- 
interacting single-particle density of states does not play a dominant role. 
At small t' , Tc is reduced by self-energy effects as discussed above and for 
intermediate values of t' the magnetic fluctuations are smaller and incom- 
mensurate so no singlet superconductivity appears. Hence at flxed U, there 
is an optimal value of t' (frustration) for superconductivity. For dx2-y2 su- 
perconductivity in under-frustrated systems (small t') Tc occurs below the 
temperature Tx where the crossover to the renormalized classical regime oc- 
curs. In other words, in under-frustrated systems at Tc the antiferromagnetic 
correlation length is much larger than the thermal de Broglie wave length 
and the renormalized classical spin fluctuations dominate. The opposite re- 
lationship between these lengths occurs for over-frustrated systems {t' larger 
than optimal) where Tc is larger than Tx and hence occurs in a regime where 
renormalized classical fluctuations do not dominate. The two temperatures, 
Tc and Tx, are comparable for optimally frustrated systems. In all cases, at 
Tc the antiferromagnetic correlation length is larger than the lattice spacing. 

Superconductivity induced by antiferromagnetic fluctuations in weak to 
intermediate coupling has also been studied by Moriya and Ueda [17] with the 
self-consistent renormalized approach that also satisfies the Mermin- Wagner 
theorem. However, in that approach there are adjustable parameters and no 
guarantee that the Pauli principle is satisfied so one cannot be certain this is 
an accurate solution to the Hubbard model. 

That there is d-wave superconductivity in the two-dimensional Hubbard 
model has by now been seen by a number of different approaches: variational 
"^^ [114,115], various Quantum Cluster approaches'^ [80-85] functional renor- 
malization group [116], and even at asymptotically small U by renormaliza- 
tion group [117]. The retardation that can be observed even tells us that spin 
fluctuations remain important for d-wave superconductivity even at strong 
coupling [77, 118, 119]. The most serious objection to the existence of d-wave 
superconductivity comes from a variational and a gaussian Quantum Monte 
Carlo approach in Ref. [120]. It could be that d-wave superconductivity in the 
two-dimensional Hubbard model is not the absolute minimum but only a local 



^^ See contribution of M. Randeria in this volume. 
^^ See contribution of D. Seneclial in this volume. 
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one. If this were the case, one could conclude that a small interaction term 
is missing in the Hubbard model to make the d-wave state the ground state. 
All other studies show that the physical properties of that state are very close 
the those of actual materials. 



1.4 More insights on the repulsive model 

The following two sections of this Chapter give a short summary of other 
results obtained with TPSC. The purpose is to show what has already been 
done, leading to the last section that contains a few open problems that could 
possibly be treated with TPSC. These two sections have some of the flavor 
of a review article, not that of a pedagogical introduction. In addition, an 
important aspect of a real review article is missing: there are very few refer- 
ences to the rest of the literature on any given topic. For anyone interested in 
pursuing some of these problems, the citation index is highly recommended. 

1.4.1 Critical behavior and phase transitions 

The self-consistent renormalized approach of Moriya-Lonzarich-Taillefer [18, 
19] was one of the first ones to treat the Hubbard model in two dimensions in a 
way that satisfies the Mermin- Wagner theorem. Other approaches exist for the 
half-filling case: Schwinger bosons [121] or constrained spin- waves [122]. The 
drawback of Moriya's approach is that it contains several fitting parameters. 
Kanamori screening, discussed in Sect. 1.2. 2, is put by hand, as is the value 
of the mode coupling constant. In addition, nothing guarantees the Pauli 
principle. In other words, Moriya's approach has much of the same physics 
as TPSC but it cannot be considered an accurate solution to the Hubbard 
model. There is also no prescription to compute the self-energy in a way that 
is consistent with double occupancy. 

More generally, the question that arises with TPSC is whether it predicts 
the correct universality class. It was shown in Ref. [38] that its results are in 
the universality class of the spherical model, namely O {N = oo) instead of 
O {N = 3) as it should be for the Hubbard model with spin-rotation invari- 
ance. This result is not surprising since the self-consistency condition on dou- 
ble occupancy found from the local moment sum rule Eq. (1.16) is very similar 
to the self-consistency condition for the spherical model. With the standard 
convention for critical exponents, one finds j/iy — 2, z — 2, and for dimen- 
sion d such that the condition 2 < d < 4 is satisfied, we find ly = 1/ {d — 2) . 
This gives in d = 3, i^ = 1,'y = 2,(3 — 1/2,1] = and 6 = 5. This should be 
compared with numerical results [123] for the 3D Heisenberg {n = 3) model, 
1^ = 0.7 and 7 = 1.4. Clearly, too close to the critical point, or too deep in the 
renormalized classical regime in d = 2, TPSC looses its accuracy. Results in 
three dimensions can be found in Refs. [55, 124]. 
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Crossover to 3d [38] 

The crossover from two- to three-dimensional critical behavior of nearly an- 
tiferromagnetic itinerant electrons was also studied in a regime where the 
intcrplane single-particle motion of electrons is quantum mechanically inco- 
herent because of thermal fluctuations. The universal renormalized classical 
crossover function from d = 2 to d = 3 for the susceptibility has been explic- 
itly computed, as well as a number of other properties such as the dependence 
of the Neel temperature on the ratio between hopping in the plane t,, and the 
hopping perpendicular to it, t±, 



1 t; 



In'^' 



(1.57) 



with Umf,c = 2/x (Qd=2,0) at the d — 2 pseudogap temperature [38]. 

Quantum critical behavior 

At r = 0, half-filling, the ground state has long-range antifcrromagnetic order. 
As one dopes, the order becomes incommensurate and eventually disappears 
at a critical point that is called "quantum critical" because it occurs at T = 
0. Such quantum critical points are common in hcavy-fermion systems for 
example. One of the surprising things about this critical point is that it affects 
the physics at surprisingly large T. 

The quantum critical behavior of TPSC in rf = 2 is in the z = 2 univer- 
sality class, like the self-consistent renormalized theory of Moriya [39]. Like 
that theory, it includes some of the logarithmic corrections found in the renor- 
malization group approach [125]. In addition, TPSC can be quantitative and 
answer the question, "How far in T does the influence of that point extend?" 
It was found [40, 126] by explicit numerical calculations away from the renor- 
malized classical regime of the d = 2 Hubbard model that logarithmic cor- 
rections are not really apparent in the range O.Olt < T < t and that the 
maximum static spin susceptibility in the (T, n)-plane obeys quantum criti- 
cal scaling. However, near the commensurate-incommensurate crossover, one 
finds obvious non-universal T and filling n dependence. Everywhere else, the 
(T, n)-depcndcnce of the non-universal scale factors is relatively weak. Strong 
deviations from scaling occur at T of order t, the degeneracy temperature. 
That high temperature limit should be contrasted with J/2 found in the 
strong coupling case [127]. In generic cases the upper limit T ~ i is well- 
above room temperature. In experiment however, the non-universality due to 
the commensurate-incommensurate crossover may make the identification of 
quantum critical scaling difficult. In addition, note that properties other than 
the maximum spin susceptibility may deviate from quantum critical scaling 
at a lower temperature [36]. 
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1.4.2 Longer range interactions 

Suppose one adds nearest-neighbor repulsion V to the Hubbard model. The 
TPSC ansatz Eq. (1.39a) can be generalized [128,129]. Then one needs to 
compute the effect of functional derivatives of the pair-correlation functions 
that appear, as in Sect. 1.2.5, in the calculation of spin and charge irreducible 
vertices. Since the Pauli principle and local spin and charge sum rules do not 
suffice, the functional derivatives are evaluated assuming particle-hole symme- 
try, which remains approximately true when the physics is dominated by states 
close to the Fermi surface. The resulting theory, called ETPSC, for Extended 
TPSC, satisfies conservation laws and the Mermin-Wagner theorem and is 
in agreement with benchmark quantum Monte Carlo results. This approach 
allows to reliably determine the crossover temperatures toward renormalized- 
classical regimes, and hence, the dominant instability as a function of U and 
V. Contrary to RPA, even the spin fluctuations are modified by the presence 
of V. Phase diagrams have been calculated. In the presence of V, charge order 
will generally compete with spin order [128, 129]. 

1.4.3 Frustration 

The ETPSC formalism outlined in the previous section is particularly im- 
portant to treat interesting problems such as that of the sodium cobaltates. 
These compounds are often modeled in an over-simplified way by the two- 
dimensional Hubbard model on the triangular lattice. To account for charge 
fluctuations, one must also include nearest-neighbor repulsion V. Even with 
this complication this is an oversimplified model. 

The density- and interaction-dependent crossover diagram for spin- and 
charge-density wave instabilities of the normal state at arbitrary wave vector 
has been computed [67]. When U dominates over V and electron filling is 
large, instabilities are mostly in the spin sector and are controlled by Fermi 
surface properties. Increasing V eventually leads to charge instabilities where 
it is mostly the wave vector dependence of the vertex that determines the 
wave vector of the instability rather than Fermi surface properties. At small 
filling, nontrivial instabilities appear only beyond the weak coupling limit. 
Charge-density wave instabilities are favored over a wide range of dopings by 
large V at wave vectors corresponding to ^3 x ^/S superlattice in real space. 
Commensurate fillings do not play a special role for this instability. Increasing 
U leads to competition with ferromagnetism. At negative values of U or V, 
neglecting superconducting fiuctuations, one finds that charge instabilities are 
favored. In general, the crossover diagram presents a rich variety of instabil- 
ities. Thermal charge-density wave fluctuations in the renormalized-classical 
regime can open a pseudogap in the single-particle spectral weight, just as 
spin or superconducting fluctuations [67]. 
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1.4.4 Thermodynamics, conserving aspects 

Conserving approaches are very popular. FLEX [15, 16] is an example. These 
approaches are attractive because they guarantee that if one evaluates the 
same physical quantity directly from the Green function or from a derivative 
of the free energy, the answer will be identical. All that is needed for a "con- 
serving" approximation is that the self-energy be generated from a Luttinger 
Ward functional [130] that enters the expression for the free energy. In addi- 
tion, conservation laws will be satisfied in transport if irreducible vertices are 
obtained from functional derivatives of the self-energy This gives a so-called 
# derivable theory Since in perturbation theory there is an infinite number 
of possible such Luttinger- Ward functionals, depending on which closed two- 
particle irreducible diagrams constructed from G and the bare interaction one 
wishes to keep, the constraint of being conserving is not a very restrictive one. 
Conserving approximations do not satisfy the Pauli principle in general and 
they sometimes give negative values of double occupancy [55] . Various other 
limitations of conserving approaches are discussed in Appendix E of Ref . [7] . 
TPSC is not obtained from a functional derivative of a Luttinger Ward 
functional. Is this a drawback? We have seen that it satisfies conservation laws 
for spin and charge at the first step. The question is whether one can find a 
unique free energy that is consistent with the one-particle Green function and 
collective modes that TPSC focuses on. This question was addressed in the 
MSc [131] and PhD thesis [126] of Sebasticn Roy. The results are summarized 
below. We conclude with an example of thermodynamic calculation in the 
context of cold atoms. 

Thermodynamic consistency 

We should really distinguish conservation laws and thermodynamic consis- 
tency. These two notions are sometimes confused, as outlined in the previous 
paragraph. We call an approach thcrmodynamically consistent when all pos- 
sible ways of computing the same thermodynamic quantity give the same 
result. 

Obtaining the self-energy from a functional derivative of the Luttinger- 
Ward functional leads to thermodynamic consistency. In TPSC there is a 
change in perspective. Instead of looking for an approximation for the free 
energy and then deducing everything else consistently, we find a single particle 
Green's function 6**^^^ or equivalently Z''^^\ as well as double occupancy 

lTr(i:(2)G(^))^C/(n,n,)(^). (1.58) 

and deduce everything else: The free energy from integration and the ir- 
reducible vertices for transport quantities from functional derivatives of 
r(2)(Sect. 1.4.5). 
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There are three ways to extract the free energy from integration: (a) cou- 
pling constant integration of double occupancy (b) integration of /u(n) and (c) 
integration of the specific heat calculated from the total energy. With the free 
energy, all thermodynamic quantities can be obtained. If we make sure that 
the three ways to compute the free energy give the same result, then there is 
thermodynamic consistency. 

With the above expression for double occupancy, Eq.(1.58), the three dif- 
ferent ways of obtaining the thermodynamic quantities are all based on the 
same object G^^^. If G*^^) were the exact solution, they would have to be 
consistent. However G^^^ is approximate. Since G^^^ satisfies all the require- 
ments for a physical Green function, it is likely to be the exact solution of 
some Hamiltonian H that is close to, but slightly different from, the Hubbard 
model. For example, H could have longer range interactions. In deriving the 
formulas for the free energy, we assume that we are working with the Hub- 
bard model. Hence, there is no guarantee that all three methods of obtaining 
F will give the same result. One can check this numerically in principle. A 
simpler test, admittedly less stringent, is to compare n (T, /x, U) obtained from 
derivatives of the three different F. 's For the nearest-neighbor hopping model 
with /3 = 10, C/ = 4 for example, the results are identical in the percent range, 
except deep in the renormalized classical regime close to half-filling where 
TPSC anyway fails [126,131]. 

The specific heat was calculated for the nearest-neighbor hopping Hubbard 
model at half-filling as a function of temperature. TPSC reproduces the peak 
observed at small temperature in QMC [132]. It is associated with the entrance 
in the renormalized classical regime. Physically, the low temperature peak is 
a remnant of the specific heat jump that would occur at finite temperature in 
mean-field theory. 

Cold atoms, entropy 

In the context of cold atoms on optical lattices, adiabatic cooling can be used 
to reach interesting low T regimes such as the pseudogap or ordered phases by 
manipulating the scattering length or the strength of the laser-induced lattice 
potential. TPSC has been used [28], and compared with QMC calculations, to 
provide isentropic curves for the two- and three-dimensional Hubbard models 
at half-filling. Since double occupancy D is extremely accurate in TPSC, the 
entropy S was computed by integrating the Maxwell relation 

i)..-(i)„,„ 

with S (T, U — Q) the known constant of integration. 

The main findings are that adiabatically turning on the interaction in d = 2 
to cool the system is not very effective. In three dimensions, adiabatic cooling 
to the antiferromagnetic phase can be achieved in such a manner, although 
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the cooling efficiency is not as tiigli as initially suggested by dynamical mean- 
field theory [133]. Adiabatic cooling by turning off the repulsion beginning at 
strong coupling is possible in certain cases. 

1.4.5 Vertex corrections and conservation laws 

Using the functional derivative formalism of Baym and Kadanoff, [58] it is 
possible to find an expression within TPSC for the optical conductivity that 
satisfies conservation laws and hence the f-suni rule. In other words, one can 
include vertex corrections. Note that the f-sum rule this times involves the 

(2) 

momentum distribution n^, obtained from the best self-energy. 

The two types of vertex corrections that are found [36] are the antifer- 
romagnetic analogs of the Aslamasov-Larkin and Maki-Thompson contribu- 
tions of superconducting fluctuations to the conductivity but, contrary to the 
latter, they include non-perturbative effects. The calculations are impossible 
unless a number of advanced numerical algorithms are used. Take the case 
with nearest-neighbor hopping only. [36] In the pseudogap regime induced 
by two-dimensional antiferromagnetic fluctuations, the effect of vertex cor- 
rections is dramatic. Without vertex corrections the resistivity increases as 
we enter the pseudogap regime. Instead, vertex corrections lead to a drop in 
resistivity, as observed in a number of high temperature superconductors. At 
high temperature, the resistivity naturally saturates at the loffe-Regel limit. 
At the quantum critical point and beyond, the resistivity displays both linear 
and quadratic temperature dependence. The disappearance of superconduc- 
tivity in the over-doped regime is correlated with the disappearance of the 
linear term in the T dependence of the resistivity [134, 135]. The relation to 
the physics of hot spots and results for other band structures (i' ^ 0) should 
appear soon. 



1.5 Attractive Hubbard model 

Working in Nambu space and following a formal procedure analogous to that 
explained in Sect. 1.2.5, one can derive TPSC for the attractive Hubbard 
model [20,57,136]. The irreducible vertex Upp in the particle-particle singlet 
channel is given by 

Upp (nj ((1 - n^)) = U (n^ (1 - n^)) (1.60) 

and is determined self-consistently at the two-particle level by the local pair 
sum-rule 

{n,n^) = {A^A) = J Vx^'^l?) cxp(-*c.„0-) (1.61) 



g 



with 
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X^^^(.)^-f%- (1.62) 

1 + ^.pXo (<?) 

and the irreducible particle-particle susceptibility 

fe 

Again the Pauli principle and a number of crucial sum rules are satisfied. So 
is the Mermin- Wagner theorem. 

In the second step of the approximation, an improved expression for the 
self-energy is obtained by using the results of the first step in an exact expres- 
sion for the self-energy, to obtain. 



Si^\k) = Un^,-U-J2Uppx'p\q)G'Lliq-k), (1.64) 



where q = (iw„,q). This is a cooperon-like formula. The required vertex 
corrections are included as required by the absence of a Migdal theorem. 
Comparisons with other approaches can be found in Rcf. [137]. 

Pseudogap from superconductivity in attractive Hubbard model 

Using the TPSC for the attractive Hubbard model, quantitative agreement 
with Monte Carlo calculations is obtained for both single-particle and two- 
particle quantities [50]. As discussed for the repulsive case in Sect. (1.2.6) one 
obtains a pseudogap in both the density of states and the single-particle spec- 
tral weight ^'^ below some characteristic temperature T* . It was even checked 
in QMC calculations that the ratio of the thermal de Broglie wavelength to 
the pairing correlation length must be larger than unity to observe the pseu- 
dogap [140]. The pseudogap, also found in Ref. [141] for example, reflects 
precursors of Bogoliubov quasiparticles that are not local pairs, contrary to 
what is often discussed in the context of the crossover from BCS to Bose 
Einstein condensation [142]. 

With increasing temperature the spectral weight fills in the pseudogap 
instead of closing it [50] . This type of behavior is obtained in high-temperature 
superconductors. The pseudogap appears earlier in the density of states than 
in the spectral function. A characteristic behavior observed at strong coupling 
appears already in TPSC at weak to intermediate coupling, namely, small 
temperature changes around T* can modify the spectral weight over frequency 
scales much larger than temperature [50]. 

Our earlier discussion about Kosterlitz-Thouless physics in Sect. (1.3.2) is 
valid in this case as well. In the attractive Hubbard model, the supercon- 
ducting transition temperature has a dome shape because at half-filling the 



^■^ For a pseudogap in the single-particle spectral weight, it is important not to 
assume a Migdal theorem, [138, 139] and include vertex corrections [7]. 
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symmetry is S0(3) so the Tc there vanishes while it is given by the finite 
Kosterhtz-Thouless Tc [143] elsewhere. The pseudogap temperature on the 
other hand decreases monotonically from half- filling where it is largest. This 
exemplifies the fact that symmetry and dimension are important to under- 
stand pseudogap physics at weak to intermediate coupling [140]. 



1.6 Open problems 

At weak to intermediate coupling, TPSC gives the best agreement with bench- 
mark Quantum Monte Carlo methods. Its strength, compared with all other 
methods, resides in a non-perturbative treatment of the Hubbard model that 
satisfies the Pauli principle and the Mermin- Wagner theorem, in addition to 
a number of other exact constraints. Also, one works in the infinite size limit 
so the effect of collective fluctuations is not limited to a small lattice like in 
dynamical mean-field theory [26] and its generalizations ^^ [144,145]. 

The main weakness of TPSC is the difficulty to extend the method be- 
yond the one-band Hubbard model. One needs to find enough sum rules to 
determine the irreducible vertices. This can be seen as a challenge and an 
opportunity for creativity. 

For example, to include nearest-neighbor Coulomb repulsion, one needs a 
way to evaluate functional derivatives of pair correlation functions to obtain 
irreducible vertices. It has been possible achieve this [129], as discussed in 
Sect. (1.4.2), but every new problem is different. As another example, take the 
case of more than one band. Then the irreducible vertices become a matrix in 
band index and one does not have enough obvious sum rules to evaluate all 
the matrix elements [146]. 

In the presence of, say, antiferromagnetism the number of irreducible ver- 
tices also multiplies and one faces the same type of challenge. To treat long- 
range ordered states with TPSC, it might be easier to start with simpler bro- 
ken symmetries such as ferromagnetism or the Pomeranchuk instability [147]. 
The interest of treating long range order is clear. For example, the renormal- 
ized classical regime of antiferromagnetic fiuctuations is presently inaccessible 
if T is much smaller than T* . Starting from the ordered state may offer and 
alternative [62]. 

The question of the interplay of disorder and interactions is a difficult but 
topical one. Far from the Anderson disorder-induced metal-insulator transi- 
tion, the impurity averaging technique [148] may prove a useful way to intro- 
duce disorder in TPSC. One may then answer the question of what happens to 
the ^ > ^th criterion for pseudogap when the mean-free path becomes shorter 
than the thermal de Broglie wave length. 

Climbing the ladder of difficult problems, the case of strong coupling [26, 
144, 145] is a real challenge. At strong coupling the self-energy is singular at 



^* See contributions of D. VoUhardt, D. Senechal, M. Potthoff and M. Jarrell in this 
volume. 
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small frequencies. In fact it diverges as l/w at half-filling. This is inconsistent 
with the starting point of TPSC where the self-energy is constant. Perhaps 
there is a way to start from the self-energy in the atomic limit inspired by 
methods that allow for multiple poles to zeroth order [12], or some other 
way [149], but it is an unsolved problem for now. 

Some problems, by contrast, appear straightforward but they can be very 
tedious. For example, in the presence of incommensurate magnetic fluctua- 
tions where singlet d-wave pairing does not occur, can triplet pairing take 
over? One can proceed along the lines of the derivation for d-wave supercon- 
ductivity [110] but with matrix source fields to generate triplet pairing in the 
Nambu formalism. The irreducible pairing vertex would again be obtained 
from functional derivatives of a matrix S'-^K 

Paring in the attractive Hubbard model is much more straightforward, as 
we have seen. It appears at the first step, without the need to generate irre- 
ducible vertices from S^^\ Nevertheless, to study the triplet channel in the 
attractive Hubbard model, one needs to introduce near-neighbor attraction V. 
That leads to the problems mentioned above in the repulsive case. As a cu- 
riosity, one could also investigate whether functional derivatives of S^^'^ in the 
attractive Hubbard model can mediate formation of order in some particle- 
hole channel. This might be a first step towards developing a method to take 
into account the different channels on the same footing in TPSC, [150] as 
is done in renormalization group approaches [151,152]. It has been found re- 
cently within the renormalization group that in quasi one-dimensional systems 
there is a strong interference between antifcrromagnetism and unconventional 
superconductivity [100]. 

Proceeding along the lines of Ref. [36] for the conductivity, it is also clearly 
possible to compute other transport quantities, such as the thermopower, but 
it is a serious computational challenge. 

There are roadblocks, but there are also opportunities for original solutions 
and breakthroughs. 
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